###########################################
# COVID
###########################################

rm (list = ls())

library(stringi)
library(Hmisc)
library(readstata13)
library(foreign)
library(lmtest)
library(sandwich)
library(stargazer)
library(msm)
library(ggplot2)
library(did)
library(lfe)
library(gridExtra)

# Load data
load("covid_data.RData")

##############
# Table A20
##############

m1 = felm(TASA_100K  ~ t_ind | REGION | 0 | REGION, data = covid, keepCX = TRUE)
m1 = tidy(m1,conf.int = T)
m1
